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In the previous paper of this series [D. P. Varn, G. S. Canright, and J. P. Crutchfield, Physical 
Review B, submitted] we detailed a procedure — e-machine spectral reconstruction — to discover and 
analyze patterns and disorder in close-packed structures as revealed in x-ray diffraction spectra. We 
argued that this computational mechanics approach is more general than the current alternative 
theory, the fault model, and that it provides a unique characterization of the disorder present. We 
demonstrated the efficacy of computational mechanics on four prototype spectra, finding that it was 
able to recover a statistical description of the underlying modular-layer stacking using e-machine 
representations. Here we use this procedure to analyze structure and disorder in four previously 
published zinc sulphide diffraction spectra. We selected zinc sulphide not only for the theoretical 
interest this material has attracted in an effort to develop an understanding of polytypism, but 
also because it displays solid-state phase transitions and experimental data is available. With 
the first spectrum we find qualitative agreement with earlier fault-model analyses, although the 
reconstructed e-machine detects structures not previously observed. In the second spectrum, the 
results cannot be expressed in terms of weak faulting and so no direct comparison between the 
fault model and the reconstructed e-machine is possible. Nonetheless, we show that the e-machine 
gives substantially better experimental agreement and a number of structural insights. In the third 
spectrum, the fault model fails completely due to the high degree of disorder present, while the 
reconstructed e-machine reproduces the experimental spectrum well. In the fourth spectrum, we 
again find good quantitative agreement with experiment but find that the e-machine has difficulty 
reproducing the shape of several Bragg-like peaks. We discuss the reasons for this. Using the e- 
machines reconstructed for each spectrum, we calculate a number of physical parameters — such as, 
stacking energies, configurational entropies, and hexagonality — and several quantities — including 
statistical complexity and excess entropy — that describe the intrinsic computational properties of 
the stacking structures. 
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I. INTRODUCTION 

In the first paperi of this two-part series we presented 
a novel technique for the discovery and description of 
planar disorder in close-packed structures (CPSs): e- 
machine spectral reconstruction theory (eMSR or "emis- 
sary"). We showed that the technique allows one to build 
the unique, minimal, and optimal model (an e-machine) 
of a material's stacking structure from diffraction spectra. 
In this sequel we demonstrate the technique using diffrac- 
tion spectra for disordered, polytypic zinc sulphide. Since 
the discovery of polytypism in mineral ZnS crystals by 
Frondel and Palache-^ in 1948, much theoretical and ex- 
perimental effort has been expended to understand this 
phenomenon. (Sec, for example, Steinberger, 4 Mardix, 5 
and Sebastian and Krishna. 6 ) ZnS is an attractive system 
to study for a number of reasons: 

(i) Simplicity of the unit cell and stacking rules. While 
many materials are known to be polytypic)^ the con- 
stituent modular layers^ (MLs) can have a complicated 



structure and complex stacking rules&2*i°« For instance, 
in ideal micas there are more than a dozen atoms in a 
unit cell, and there are six ways two MLs can be stacked. 
Kaolins and cronstedtites present even more complex- 
ity^ In contrast, ZnS is simple in the extreme: its basis 
is composed of but two atoms — a zinc and a sulphur. 6 
They are arranged in a double close-packed hexagonal 
net, with one species displaced relative to the other by a 
quarter body diagonal (as measured by the conventional 
unit cubic cell) along the stacking directioniAi We take 
this double close-packed layer to be a ML£ The stacking 
of MLs proceeds as for all CPSs; namely, there are three 
absolute orientations each ML can occupy — call them A, 
B, and C — with the familiar stacking constraint that no 
two adjacent MLs have the same orientation^ 

(ii) Complex polytypism. ZnS is one of the most poly- 
typic substances known with over 185 identified crys- 
talline structures Of these, only about a dozen fairly 
short-period polytypes (up to 21 MLs) are common in 
mineral ZnS, with the remainder found in synthetic crys- 
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tals. Some of the crystal structures have repeat distances 
that extend over 100 MLs. Also, many structures show 
considerable planar disorder. The wide diversity of struc- 
tural complexity remains one of the central mysteries of 
polytypism. 

(iii) Solid-state transformations. It is believed that there 
are only two stable phases of ZnS, the low-temperature 
modification being the /3-ZnS or sphalerite (3G— ) and 
the high-temperature modification wurtzite (2H) or a- 
ZnS£ The former transforms enantiotropically into the 
latter at 1024 C. The plethora of structures suggests that 
most of them are not in equilibrium but rather structures 
that are trapped in a local minimum of the free energy 
and lack the necessary activation energy to explore all 
of configuration space. It is possible to observe these 
structures by annealing and then arresting the transfor- 
mation upon quenching. One can then study the various 
intermediate stages of the transformation. 

(iv) Availability. Polytypes of ZnS, both ordered and dis- 
ordered, are easily manufactured in the laboratory by a 
variety of methods^ One of the more common is growth 
from the vapor phase above temperatures of about 1100 
C. Crystals can also be grown from melt at high pres- 
sures, by use of chemical transport and hydrothermally. 
The distribution of polytypes observed depends on the 
method used. 

Nearly a dozen theories have been proposed to explain 
polytypism, among them being the ANNNI modeljiii^ 
Jagodzinski's disorder theory, 15 and Frank's screw dis- 
location theory. 16 (For a complete discussion, see for ex- 
ample, Verma and Krishna^ Trigunayat^ and Sebastian 
and Krishna^) We will have little to say here about the 
mechanisms that produce various polytype structures. 
Instead, our focus will be on describing the disordered 
structures so commonly seen. We feel that an adequate 
description of the disordered structures — which so far has 
been lacking — is warranted before one can evaluate mod- 
els that explain the formation of disorder and structures 
and, especially, the solid-state phase transitions that lead 
to them. 

Previous descriptions of planar disorder in single 
crystals of ZnS fall into one of two categories: the 
fault model (FM)ii& and Jagodzinski's disorder model 
(DM). 19 i 20 i 21 i 22 Applications of the FM include Roth's 2 ^ 
study of faulting induced in hexagonal crystals grown 
from the vapor phase upon annealing. Roth extracted 
correlation information from the diffraction spectra by 
Fourier analysis and then derived analytical expressions 
relating how correlation functions decayed with both in- 
creasing separation between MLs and as a function of 
the fault probability. He considered both randomly dis- 
tributed growth and deformation faults and found that 
for weakly disordered specimens deformation faulting 
gave the best agreement with experiment. 

Significant applications of the FM to planar disorder 
in ZnS have been carried out by Sebastian, Krishna, and 



coworkers. They studied the 2H-to-3C solid-state trans- 
formation in vapor-grown ZnS crystals after annealing 
between temperatures of 300 to 650 C~ By analyzing 
and comparing the profiles of the integer-Z reflections 25 
to those of the half-integer-i reflections for weakly faulted 
crystals, they found that the disorder was largely due to 
the random insertion of deformation faults. They at- 
tributed slight discrepancies between the observed and 
calculated profiles to the so-called nonrandom insertion 
of faults. Sebastian and Krishna 2 — later studied the disor- 
dered stacking in 3C crystals grown from the vapor phase, 
as well as those obtained from annealing 2H crystals. 
They found that the structure of both the as-grown and 
annealed crystals was best explained as randomly dis- 
tributed twin faults in the 3C structure. They concluded 
that the 2H-to-3C transformation in ZnS proceeded by 
the nonrandom nucleation of deformation faults occur- 
ring preferentially at two ML separations. 

To better understand the nature of the nonrandom 
insertion of deformation faults in the 2H structure, Se- 
bastian and Krishna-2i2& introduced a three-parameter 
model that assigned separate probabilities to the random 
insertion of deformation faults, as well as deformation 
faulting at two and three ML separations. They derived 
analytical expressions for the diffraction spectra in terms 
of these parameters and concluded that both the 2H-to- 
3C and the 2H-to-6H transformations proceeded via the 
nonrandom nucleation of deformation faults. Their anal- 
ysis showed that these transformations occurred simul- 
taneously in different regions of the same crystal. They 
attributed this to variations in the stoichiometry. Sebas- 
tian 2 ^ gave a similar treatment that came to the same 
conclusions. With the exception of Roth, all of these 
analyses depended on carefully characterizing the change 
in Bragg peaks as one introduces a small amount of dis- 
order. We have previously given a criticism of the FM 
approach elsewhereii^ 

Jagodzinski's DM is a two-parameter model that as- 
sumes two thermodynamically stable phases in CPSs: the 
2H and 3C. One therefore finds two kinds of fault (and 
here we mean structure and not mechanism): namely, 
cubic faults in the 2H and hexagonal faults in the 3C. By 
choosing appropriate values of the two model parameters 
one can also model 4H structures. Within this frame- 
work, an analytical expression for the diffracted inten- 
sity is derived that depends on the model parameters 
in a complicated manner. Nonetheless, one can select 
model parameters that give the best agreement with ex- 
periment. 

Muller— used this method to analyze faulted ZnS 
diffraction spectra and found that while he was able to 
obtain reasonable agreement between theory and exper- 
iment for a few spectra, for many he was not. Singer 31 
re-examined this approach and concluded that the DM 
applies when faulting is random, but when the faulting 
is nonrandom, as many ZnS specimens are suspected to 
be, the model fails. However, Frey et al22k studied the 
3C-to-2H transformation in single crystals of ZnS using 
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the DM and were able to obtain excellent agreement be- 
tween theory and experiment. They fitted the experi- 
mental diffraction spectra to the DM's analytical one. 
Due to the complicated nature of the expression, how- 
ever, they treated eight constants that depend on the 
two model parameters as independent. From these eight 
fitted parameters they were able to find the two model 
parameters that best fit each spectra. 

One cannot help but raise questions concerning the 
mathematical rigor used to find the model parameters 
in this way. We show elsewhere 3 ^ that the description 
of the stacking disorder as given by the DM is a special 
constrained case of the r — 2 computational mechanics 
approach. As in the latter, in the DM there is no assump- 
tion of weak faulting and one does use diffuse scattering 
to build the model. Since the spectra we analyze have 
not been previously treated using the DM and it is a spe- 
cial case of our own, we do not discuss the DM further 
here, but treat it in detail elsewhere* 3 ^. 

A third possible method of discovering structural infor- 
mation about disordered solids from diffraction spectra 
employs a reverse Monte Carlo (RMC) techniquei^24 In 
this method, one typically searches for a configuration of 
constituent atoms such that a signal — e.g., the diffrac- 
tion spectrum — estimated from the candidate structure 
most closely matches the experimental signal. This tech- 
nique can be applied for the case of disorder in three 
dimensions. One drawback, however, is that candidate 
structures are often found that are physically implausi- 
ble. One needs to impose assumptions to eliminate these. 
To our knowledge, this technique has not been applied to 
polytypism, and we do not consider it here. 

In this work, we apply computational mechanics 35 - 36-37 
to discover and describe disordered stacking sequences in 
four previously published ZnS diffraction spectra. We 
Fourier analyze each spectrum to find correlation infor- 
mation between MLs and then calculate the probability 
distribution of stacking sequences. From the latter, we 
reconstruct the e-machine that gives the stochastic pro- 
cess for the ML stacking and compare it to previous FM 
analyses. From the reconstructed e-machine, we calculate 
the stacking entropy per layer, average stacking-fault en- 
ergy per Zn-S pair, memory length, hexagonality, and 
generalized period^ We find that the diffraction spectra 
of the four samples is well described using the computa- 
tional mechanics approach. 

We note that our primary purpose in the following is 
expository; that is, we wish to demonstrate the efficacy 
of eMSR on real data. Since we use diffraction spectra 
from older studies) 6 - the analyses given here are less than 
ideal. Specifically, we digitized data from the published 
spectra and found that there was significant systematic 
error in each spectrum. Additionally, the experimental 
data was not reported with error bars. Despite these pos- 
sible shortcomings, eMSR allows us to offer more com- 
prehensive analyses of the spectra than given by previous 
workers. 



Our development is organized as follows. In fjn]we out- 
line our approach, including a brief discussion of the ex- 
perimental methods and our analysis. In 31111 we give the 
results of e-machine reconstruction for four experimental 
ZnS diffraction spectra and contrast this to the FM ap- 
proach when possible. In fllVI we calculate the stacking 
energies per Zn-S pair and the hexagonality for the vari- 
ous structures from our reconstructed e-machines. In §V\ 
we give our conclusions and propose some directions for 
future theoretical and experimental work. 



II. METHODS 

The four diffraction spectra we analyze come from Se- 
bastian and Krishna^ and are labeled SKXXX by the 
page (XXX) on which they appear in that source. These 
data were collected in the 1980's and since they no longer 
exist in numerical form^ we digitized them from the 
diffractograms given in the Sebastian and Krishna pub- 
lication^ In this section, we give a brief synopsis of the 
experimental procedure, discuss the assumptions made 
to analyze the data, and list the corrections we apply to 
the experimental spectra. 



A. Experimental Details 

The experimental procedure is given in more detail 
elsewhere Briefly, the crystals were grown from 
the vapor phase at a temperature in excess of 1100 C in 
the presence of H2S gas. Each crystal was needle-shaped, 
approximately 0.1 to 0.4 mm in diameter and 1 to 2 mm 
in length. Two of the four crystals were further annealed 
for one hour at 300 and 500 C. These investigations were 
performed to better understand the fault structures they 
contain, as well as study the solid-state transformations 
that ZnS crystals undergo. 

The intensity along the 10.Z reciprocal lattice row was 
recorded using a four-circle single-crystal diffractometer 
for each specimen in steps of approximately = 0.005. 
(Our definition of I differs slightly from that of Sebastian 
and Krishna^ so the ^-increment we report also differs.) 
The crystal and the counter were held stationary while 
the crystal was illuminated with WLoK a radiation. The 
sharp reflections along the h — k = mod 3 rows were 
used to orient each crystal. The divergence of the inci- 
dent beam was adjusted to cover the mosaic spread for 
each crystal. The experimental diffraction spectrum is 
reported as the total number of counts versus I. The 
crystals were examined under a vickers projection mi- 
croscope and did not show signs of kinking or shearing, 
even after annealing. They did show parallel striations 
or stripes perpendicular to the stacking direction. 
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B. Assumptions 

To make the analysis tractable, we employ the follow- 
ing assumptions common in the analysis of planar disor- 
der in ZnS: 

(i) Each ML is perfect and free of distortions and defects. 
We assume that each ML is identical and the MLs them- 
selves are undefected. That is, each ML is crystalline 
in the strict sense, with no point defects, impurities, 
or distortions in the two-dimensional lattice structure. 
This clearly precludes the possibility of screw disloca- 
tions, which are known to play a role in the polytypism 
of some ZnS crystals Since each of the crystals we 
analyze was examined under a vickers projection micro- 
scope and no such dislocations were seen, and the crystals 
retained their shape after annealing, this seems reason- 
able. (It is known 6 , that during solid-state transforma- 
tions of specimens of ZnS with an axial screw dislocation 
the specimen will exhibit " kinking" with a characteristic 
angle of 19°28'.) 

(ii) The spacing between MLs is independent of the local 
stacking arrangement. There is known to be some slight 
dependence of the inter-ML spacing depending on the 
polytype&2S For the 2H structure in ZnS, the inter-ML 
spacing is measured to be 3.117 A. For the 3C structure, 
the cubic cell dimension is a c — 5.412 A which gives a 
corresponding inter-ML spacing of 3.125 A. Therefore, to 
an excellent approximation, this spacing is independent 
of the local stacking environment in ZnS. 

(iii) The scattering power of each ML is the same. We 
assume that each ML diffracts x-rays with the same in- 
tensity. There is no reason to believe that this is not so, 
unless absorption effects are important or the geometry 
of the crystal is such that each ML does not have the 
same cross-sectional area. 

(iv) The stacking faults extend over the entire fault plane. 
Examination under microscope indicates that this is gen- 
erally true. However, Akizuki 40 found evidence that the 
faults do not extend completely over the faulted MLs by 
examining a partially transformed ZnS crystal under an 
electron microscope. 

(v) We assume that the "stacking process" is stationary. 
We simply mean that the faults are uniformly distributed 
though out the crystal. Put another way, we assume 
that probability of finding a particular stacking sequence 
is independent of its location in the crystal. This does 
not, however, preclude regions of crystal structure in- 
terspersed between regions of disorder. It simply means 
that the statistics of the stacking does not change as one 
moves from one end of the crystal to the other. 

Notably absent from this list are any assumptions 
about the crystal structure present (if any) and how the 
sample might deviate from that structure. In contrast 
to the FM, we invoke no a priori structural assumptions 
concerning the stacking sequence. 



C. Corrections to the Experimental Spectra 

We corrected each spectrum for the following effects: 

(i) The atomic scattering factor. This correction ac- 
counts for the spatial distribution of electrons, as well 
as for the wavelength of the incident radiation and angle 
of reflection. Calculations of these effects are given in 
standard tables 41 and we employ them in our work. 

(ii) The structure factor. The structure factor— accounts 
for the two-atom basis in ZnS. 

(iii) Anomalous scattering factors. Also called dispersion 
factors, the anomalous scattering factors correct for the 
binding energy of the electrons and the phase shifts4i 
For our case, we find these to be small, but have included 
them nonetheless. 

(iv) Polarization factor. We use the standard correction 
factor for unpolarized radiation4£ 

Factors we did not correct for include the following: 

(i) Thermal factors. At room temperature, this effect is 
small for ZnSi£ 

(ii) Absorption factor. For the geometry of the ZnS crys- 
tals we analyze, the linear mass coefficient^ 4 ^ is much 
larger than the thickness, therefore we ignore it. 

(iii) Instrument resolution. This is not reported with the 
experimental data, so we do not deconvolve the spectrum. 



III. ANALYSIS 

We now give the structural analysis of four experi- 
mental diffraction spectra taken from Sebastian and Kr- 
ishna. 6 We apply eMSR to each to build a model that 
describes the stacking process. From our model, we cal- 
culate various measures of intrinsic computation for each 
spectra. We compare our results with that obtained by 
previous researchers using the FM. Since the experimen- 
tal spectra are not reported with error bars, we are unable 
to set an error threshold T as required in eMSR. Instead, 
we found that each increase in the memory length r con- 
tinues to give a better description of each spectral We 
perform e-machine reconstruction up to r — 3 for each 
spectra. To find the CFs from the e-machines, we take 
a sample length 400 000, as generated by the e-machine. 
The diffraction spectra are calculated using 10 000 MLs. 
The experimental spectra are normalized to unity over 
the Z-interval used for reconstruction, as are the spec- 
tra calculated from each e-machine. For the spectra cal- 
culated from the FM, we set the overall scale to best 
describe the Bragg peaks as shown in Sebastian and Kr- 
ishna. 6 We also calculate the profile 72.-fact.or 1 to evaluate 
the agreement between experiment and theory for each 
spectrum. 
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FIG. 1: Diffraction spectrum along the 10. 1 row from SK134 
corrected for atomic scattering factors, the structure factor, 
dispersion factors, and the polarization factorp— We see 
that the spectrum is not periodic in I, as it should be. This 
indicates that there exist significant errors in the data. 



FIG. 2: Experimental and theoretical figures-of-merit — /3 and 
7 — as a function of lo for diffraction spectrum SK134. We 
define lo as the point at which integration over a unit ^-interval 
is initiated. We see that lo ~ 0.04 gives the best agreement 
with the theoretical values. 



A. SK134 

The corrected diffraction spectrum along the 10. 1 row 
for an as-grown 2H ZnS crystal is shown in Fig. 2] One 
immediately notices that the spectrum is not periodic in 
I, as it should be, but instead suffers from variations in 
the intensity. The peaks at I = 0.0 and I — 1.0 are of sim- 
ilar intensity, but the peaks at I — —0.5, 0.5 and 1.5 seem 
to suffer from a steady decline in intensity. We can there- 
fore be confident that this spectrum contains substantial 
systematic error as reported by the experimentalists. 

This difference in diffracted intensity between peaks 
results from the finite thickness of the Ewald sphere due 
to the divergence of the incident beamm^ A suitable 
choice of geometry can minimize, but not eliminate these 
effects, such that one finds only a gradual variation in / 
with Al. Since analysis by the FM depends only on the 
change of the shape and the position of the Bragg-like 
peaks, such a slow variation of the diffracted intensity 
with I will not affect the conclusions drawn from an FM 
analysis. It is possible to correct these effects but 
this has not been done in the literature. 

Our analysis depends on selecting an appropriate l- 
interval where variations in diffracted intensity due to 
these experimental effects are minimized. It is impor- 
tant, then, to select an interval that is relatively error- 
free. As discussed in Part Iji there are two criteria, called 
figures-of-merit and denoted 7 and (3, one can use for 
this. It can be shown that in an error-free spectrum the 
parameters must be equal to constant values —1/2 and 
1, respectively, for any unit ^-interval, regardless of the 
amount of planar disorder present. The extent that (3 
and 7 differ from their theoretical values over a given l- 
interval measures how well the diffraction spectrum over 
the interval can be represented by a physical stacking of 
MLs. It makes sense, then, to choose an Z-interval for e- 
machine reconstruction such that the theoretical values of 




FIG. 3: The r = 3 e-machine reconstructed for SK134. The 
recurrent causal states {So ■ ■ ■ S7} are labeled (in decimal no- 
tation) by the last three (binary) spins seen and the asymp- 
totic state probabilities are given for each state. Edge la- 
bels s\p indicate a transition on spin s with probability p. 
The large asymptotic state probabilities for the S2 and <Ss 
states, as well as the large state-transition probabilities be- 
tween them, show that this is predominately a 2H crystal 
with some faulting. 



the figures-of-merit are most closely realized. This does 
not, of course, guarantee that the interval is error-free. 
Glancing at Fig. [21 shows that 7 = —0.51 and [3 = 0.95 
over the interval I € [0.04, 1.04]. 

We perform eMSRi and find that the smallest-r e- 
machine that gives reasonable agreement between the 
measured and e-machine spectra has a memory length of 
r = 3. The reconstructed e-machine is shown in Fig. [3] 
The large asymptotic state probabilities for the S2 and 1S5 
causal states (CSs), as well as the large inter-state transi- 
tion probabilities between them, indicate this is predom- 
inantly a 2H crystal. More specifically, the probability 
of seeing sequences 1010 and 0101, corresponding to the 
2H cycle, have a combined total weight of about 64%. 



6 



Structure 


e-Machine Fault Model 


2H 


64% 


83% 


3C 


8% 


0% 


Deformation fault 


16% 


17% 


Growth fault 


6% 


0% 


Layer-displacement fault 


6% 


0% 



TABLE I: Structural decomposition of SK134 according to 
the reconstructed e-machine of Fig. |3] and according to the 
fault model analysis of Sebastian and Krishnai— The latter is 
valid only under the assumption of weak faulting. 



The remaining probability is distributed among the other 
fourteen length-4 sequences. It is tempting to interpret 
the remaining structure in terms of faults and, indeed, it 
seems we can. 

Let us treat the transitions s = from 1S4 and s = 1 
from So as though they are missing for the purposes of a 
fault analysis. This implies that the sequences 0001 and 
1000 are disallowed. Of course, this cannot be exactly 
true, as the CS So would then be isolated from the rest 
of the e-machine. In this case, we would say that the 
e-machine is not strongly connected and, as such, cannot 
represent a physical stacking of MLs. However, the com- 
bined probability weight of these two sequences is < 1%, 
so neglecting them gives only a small error in our intu- 
itive understanding of the faulting structure. 

Then on the left half of the e-machine there is struc- 
ture associated with a 2H deformation fault [iSa^^^s] 
with probability weight 0.16 = P(1011) + P(0111) + 
P(1110) +P(1101). We can interpret the causal-state cy- 
cle (CSC) [<S3<S6<S4<Si<S2iS5] as a layer-displacement fault 
and see that it is assigned a probability weight of 0.06. 
The right portion contains the CSC [S41SHS2] with proba- 
bility weight 0.06, which is associated with growth faults. 
The CSCs [So] and [S7], identified as 3C structure, have 
a combined weight of 0.08. 

Given these observations, a possible interpretation sug- 
gested by the e-machine is that SK134 has crystal struc- 
tures and faults in the proportions given in Tabled The 
decomposition there is sensible since there is an under- 
lying crystal structure present and the smaller, faulting 
paths are not too large. As we will see, these need not 
always be the case. Sebastian and Krishna^ have ana- 
lyzed this diffraction spectrum using the FM and found 
that approximately one in every twenty MLs is deforma- 
tion faulted, so they described the stacking structure as 
a faulted 2H crystal with 5% random deformation fault- 
ing. This is equivalent to assigning CSCs responsible for 
deformation faulting a total probability weight of 0.17. 
We compare the structure analyses of the two models in 
Tabled 

We see that both analyses agree that the dominant 
structure is 2H, though the e-machine attributes less 
of the crystal structure to this "parent" phase. Simi- 
larly, both find that structures associated with deforma- 
tion faulting are important and assign them almost equal 



weights. 

They differ, however, in that the e-machine finds addi- 
tional faulting structures (growth and layer-displacement 
faulting), as well as 3C crystal structure. Since this 
crystal, if annealed at sufficient temperatures for a long 
enough time, will transform into a twinned 3C structure, 
the latter is easily understood as nascent structure in 
that process. Finding the presence of weak 3C structure 
is not unreasonable since there is some slight enhance- 
ment of the diffracted intensity at I w 1/3 and I w 2/3. 
The other faulting structures seen are less easily under- 
stood. Growth faults, so-named because they primarily 
are formed during the growth of the 2H crystal, are not 
expected to play an important role in solid-state transfor- 
mations of ZnS, 23 Their presence here may result from 
the initial growth of the crystal. The small amount of 
layer-displacement structure could be seen as two adja- 
cent, yet oppositely oriented deformation faults. That is, 
a deformation fault in a 2H structure is simply a spin flip 
in the Hagg notation^ so that a sequence . . . 010101 . . . 
transforms to . . . 01L101 ... as a result of one deforma- 
tion and then to ...011001... upon another resulting 
in a layer-displacement fault . . . 011001 .... This might 
imply some coordination between faults. Or, the mecha- 
nism of layer-displacement faulting may play some minor 
role in the solid-state transformation. 

However, one cannot disambiguate these from the 
available spectra and the reconstructed e-machine. The 
e-machine only provides information about the structure. 
We must look outside the e-machine to formulate an un- 
derstanding of how the polytype came to be stacked in 
this way. This is the critical difference between faulting 
mechanism and faulting structure. In the former, a phys- 
ical process is responsible for causing the MLs to shift or 
deviate from a perfect crystal. In the latter, in the limit 
of weak faulting, the physical process leads to a given 
(statistical) structure. In this limit it may be possible 
to postulate with some certainty that the mechanism re- 
sulted in the observed structure. For more heavily faulted 
crystals, however, such an identification of structure with 
mechanism is dubious. Other techniques, such as numer- 
ical simulationa 47 i 48 i 49 i 50 i 51 i 52 i 53 i 54 or analysis of a series 
of crystals in various stages of the transformation, are 
necessary to unambiguously determine the mechanism. 

Returning to the analysis of SK134, Fig. 0] compares 
the CFs obtained from the experimental diffraction spec- 
trum, those obtained from the e-machine and those from 
the FM. There is reasonable agreement between the ex- 
perimental and e-machine-predicted CFs. For small n, 
however, the FM overestimates the amplitude in the os- 
cillations in Q s (n). The experimental diffraction spec- 
trum is compared to that calculated from the FM and 
e-machine in Fig. [3] Both models give good agreement 
near the Bragg peaks at I — 0.5 and I — 1.0, with per- 
haps the FM performing a little better at I = 1.0. The 
diffuse scattering near the shoulders of the I = 0.5 peak 
are better represented by the e-machine. We calculate 
the profile 7?.-factor between experiment and the FM to 



7 




2 4 6 8 10 12 14 16 18 20 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0.1 0.2 

Jl I 



FIG. 4: Q 3 (n) versus n for experimental spectrum SK134 
(open squares connected by solid line) , the fault model (open 
circles), and the r — 3 e- machine (solid squares). The Q a (n) 
are defined only for integer values of n, but lines are drawn 
connecting adjacent points as an aid to the eye. We see good 
agreement up to n ~ 15, after which the r = 3 approximate 
correlation functions die too quickly to their asymptotic value 
of 1/3. 
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FIG. 5: Comparison of the experimental diffraction spectrum 
SK134 along the 10./ row (triangles) for a disordered ZnS sin- 
gle crystal with a spectra estimated from the FM with 5% 
deformation faulting (dashed line) and r = 3 e-machine (solid 
line). The vertical scale in the inset is logarithmic intensity. 
For clarity, we report only a few representative data points 
from the experimental diffraction spectrum. In our analysis, 
however, we used the much finer mesh reported in the experi- 
mental data. We find that the eMSR gives 1Z ( m = 20%, while 
the FM gives Kfm = 33%. 



be li-FM = 33% and between experiment and eMSR to 
be TZ eM = 20%. 

From the e-machine we calculate the three length pa- 
rameters to be r = 3 ML, V = 4.8 ML, and A c = 9.5±0.5 
ML. The three measures of intrinsic computation are 
found to be = 0.50 bits/ML, C p = 2.3 bits, and 
E = 0.75 bits. We return to discuss and compare these 
for all of the samples in a later section. 



FIG. 6: Comparison of the experimental diffraction spectrum 
SK135 along the 10. 1 row (triangles) for a disordered 3C single 
crystal with the diffraction spectra calculated from the FM 
with 12% twinned faulting (dashed line) and the e-machine 
(solid line). We find that the eMSR gives Km = 13%, while 
the FM gives ft FM = 33%. 



B. SK135 

The next sample we examine is a twinned 3C crystal 
obtained by annealing a 2H crystal at 500 C for 1 h. The 
diffraction spectrum for this crystal is given in Fig. El We 
find that figures-of-merit are closest to their theoretical 
values over the interval I E [—0.80, 0.20] with values (3 = 
—0.50 and 7 = 0.93. The smallest-r e-machine that gives 
reasonable agreement with experiment was found at r = 
3 and has a profile 7\L-factor of TZ 6 m — 13%. The resulting 
e-machine is shown in Fig. [7] Based on the presence of 
asymmetrically broadened peaks and the absence of peak 
shifts, a FM analysis^ finds this sample to be a twinned 
3C crystal with 12% twinned faulting. The profile 1Z- 
factor between experiment and the FM is found to be 
Ti-FM — 33%. 

The large CS probabilities associated with CSs So and 
1S7, as well as their large self- loop transition probabil- 
ities, suggest that this is a twinned 3C crystal. Wc 
also note that the transitions corresponding to antiferro- 
magnetic paths (0101 and 1010) have a relatively small 
combined weight of only about 4%. In fact, the prob- 
ability weight for the 0101 path is zero. (The transi- 
tion from S2 to 1S5 is missing.) This indicates that 2H 
structure has largely been eliminated. In addition to 
the 0101 path, the 1001 and 0010 paths are also miss- 
ing. This implies that twinned faulting is important, but 
also the remnant of the 1010 path has some role. In- 
stead of a simple twinned fault [<S7<S6<5>4<Sb] gi ym S the 
sequence . . . 1111 10000 . . ., where the vertical line indi- 
cates the fault plane, the path [S71S655S2S5S0] giving 
the sequence . . . 1111|01000 . . . has approximately twice 
as much probability weight associated with it. In the 
e-machine's right portion twinned faulting [<So<Si«S3<S7] is 
largely responsible for the (0)* 3C cycle converting to the 
(1)* 3C cycle and we also observe that double deforma- 
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FIG. 7: The reconstructed r = 3 e-machine for SK135. The 
strong self-loop transition probabilities between causal states 
So and 57 as well as their large asymptotic state probabilities 
indicate that the . . . 0000 . . . and ... 1111 .. . structures are 
important. This, then, is a twinned 3C crystal. 



FIG. 8: Comparison of the Q s (n) versus n for experimen- 
tal spectrum SK135 (open squares), the e-machine (solid 
squares), and the FM (open circles). The e-machine gives 
good agreement with experiment, while the FM overestimates 
the oscillation magnitude. 



tion faulting [<So<5>i<S3<S6<S4] pl a Y s a r °l e - 

It is interesting to mention that, while a ML of ZnS has 
spin-inversion symmetry^ and, thus, the one-dimensional 
Hamiltonian describing the energetics of the stacking is 
also spin invariant, in general e-machines need not be 
spin-inversion invariant. [Note that the e-machine in 
Fig.UJis not spin-inversion invariant.] There is, of course, 
no reason why we should demand spin-inversion invari- 
ance. After all, then one could never have a crystal of 
purely one 3C structure or the other. However, since 
this crystal was initially in the 2H structure — which is 
spin-inversion invariant — it is curious that this is not 
preserved as the crystal is annealed. That is, there is 
no reason to expect that faulting should occur preferen- 
tially with one chirality. Notably, the FM always assumes 
spin-inversion symmetry. 

Examining the Q s (n) estimated from experiment with 
those found from the e-machine in Fig. EI we find reason- 
able agreement up to n ~ 20. The Q s (n) found from the 
FM generally overstate the magnitude of the oscillations 
in the CFs. 

We can further examine the diffraction spectra. In 
Fig. El the diffraction spectrum found from the FM and 
the e-machine are compared with experiment. The e- 
machine gives a good fit, except perhaps at a shoulder in 
the experimental spectrum at I = —0.6 and the small rise 
at I = —0.16. The latter might be understood as a minor 
competition between the 3C and 6H CSCs that is not 
being well modeled at r = 3. Comparison of the diffrac- 
tion spectrum from the FM with that from experiment 
reveals good agreement with the peak at I = —0.33 and 
poor agreement with the one at I = —0.67. This is not 
surprising as the FM did not use the peak at I = —0.67 
to find the faulting structure. Likewise, the diffuse scat- 
tering between peaks is not at all well represented by the 
FM. Additionally, the small rise in diffuse scattering at 
/ = —0.16 is likewise absent in the FM diffraction spec- 
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FIG. 9: Comparison of the experimental diffraction spectrum 
SK137 along the 10.1 row (triangles) for a disordered 3C single 
crystal with the diffraction spectra calculated from the FM 
with 6.8% twinned faulting (dashed line) and the e-machine 
(solid line) . The profile 1Z- factor between experiment and the 
e-machine calculated diffraction pattern is lZ t M = 17%. The 
FM gives considerably worse agreement, with a calculated 
profile 7?.-factor of TZfm = 58% between it and experiment. 



trum. 

From the e-machine we calculate the three length pa- 
rameters to be r = 3 ML, V = 5.6 ML, and A c = 4.4±0.7 
ML. The three measures of intrinsic computation are 
found to be = 0.59 bits/ML, C M = 2.5 bits, and 
E = 0.71 bits. 



C. SK137 

The third experimental spectrum we analyze comes 
from an as-grown disordered, twinned 3C crystal. The 
diffraction spectrum for this crystal along the 10.1 row 
is shown in Fig. |5J The figures-of-merit are closest to 
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their theoretical values over the interval I G [—0.8, 0.2] 
with values of 7 = —0.49 and [3 = 0.98. We performed e- 
machine reconstruction up to r = 3 and found that this 
produces reasonable agreement with experiment giving 
a profile 7^-factor of H e M — 17%. The r = 3 approxi- 
mate e-machine in shown in Fig. A FM analysis finds 
SK137 to be a twinned 3C crystal with 6.8% twinned 
faulting**. The FM-calculated diffraction spectrum has a 
profile 7^-factor of lZ e M = 58% when compared with ex- 
periment. 

A comparison of the CFs from experiment, the FM, 
and the e-machine is shown in Fig. 1111 For smaller n, 
the e-machine gives good agreement with experiment, al- 
though the error increases at larger n. As shown else- 
where^ the experimental CFs maintain small, but per- 
sistent oscillations about their asymptotic value of 1/3 
up to n w 40, while the CFs derived from the e-machine 
effectively reach this asymptotic value at n w 25. This 
leads us to speculate that there is some structure in the 
stacking process that the e-machine is missing. We ex- 
pect that reconstruction at r — 4 will prove interesting 
here. This has not yet been completed. 

The FM fares markedly worse. It substantially overes- 
timates the strength of the oscillations in the CFs for all 
n. 

A comparison of the diffraction spectra for experiment, 
the FM, and the e-machine is given in Fig. [5] The e- 
machine gives reasonable agreement everywhere except 
around the Bragg peaks at I = —0.33 and I = —0.67. 
Here the e-machine gives a value for the peak intensity 
15% and 35% lower, respectively, than experiment. The 
FM does much better at the Bragg peaks, as one might 
expect. The diffuse scattering between the peaks, and 
especially the broad-band rise in intensity near I = —0.5, 
are simply missing in the FM diffraction spectrum, how- 
ever. The e-machine fit in this region is substantially bet- 
ter, picking up a number of important spectral features, 
such as broadband components and broadened peaks. 

What does the e-machine imply about the stacking 
process? All CSs and allowed transitions are present ex- 
cept for the transition between S\ and 53. This absent 
transition implies that the 0011 stacking sequence is not 
present in SK137. This, then, means that the 000111 
sequence, and hence the CSC [<S7<S6<S4<So<Si«S3] associated 
with the 6H structure^ is also absent. Therefore, in this 
twinned 3C crystal there is no 6H structure. This is sur- 
prising, since many ZnS spectra show enhancement about 
the 6H positions during solid-state phase transitions from 
2H to 3C. In Fig.[!5]there is (arguably) a slight increase in 
diffracted intensity at I = —0.16 and I = 0.16. So the ab- 
sence of the 6H structure does seem echoed in the exper- 
imental spectrum. There is, however, a large broadband 
increase in intensity about I = —0.5 and a much smaller 
increase about I = 0.0. Reflections at these I are usually 
associated with 2H structure, with the half-integer peaks 
carrying three times the intensity of the integer peaks. 
The e-machine does show that the CSC [<S2<Ss] associ- 
ated with the 2H structure is present. The frequency of 




FIG. 10: The reconstructed r = 3 e-machine for SK137. The 
strong self-loop transition probabilities between causal states 
So and S7, as well as their large asymptotic state probabilities, 
suggest that the . . . 0000 . . . and ... 1111 .. . structures are im- 
portant. Notice that, unlike the e-machine for SK135, the 
CSC [iS2<Ss] is present, suggesting that associated 2H struc- 
ture is present. The absence of the transition between CSs 
<Si and <S3 implies that the 0011 sequence, and therefore the 
CSC associated with the 6H structure, is not present. 
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FIG. 11: Comparison of the Q 3 (n) versus n for experimental 
spectrum SK137 (open squares), the r — 3 approximate e- 
machine (solid squares), and the FM (open circles). 



occurrence of the 0101 and 1010 stacking sequences to- 
gether make up about 12% of the total probability weight 
on the e-machine. Even though Pcsc {[S2S5]) <C 1, it is 
not unreasonable to suggest that 2H structure is present. 

As with the other spectra, we can calculate 
several characteristic length and information- and 
computational-theoretic quantities from the e-machine. 
We find the three length parameters to be r — 3 ML, 
V = 6.7 ML, and A c = 12 ± 3 ML. The three measures 
of intrinsic computation are found to be = 0.65 bits, 
Cf, = 2.7 bits, and E = 0.79 bits. 
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FIG. 12: A comparison of the experimental diffraction spec- 
trum SK229 (dashed line) and that calculated from the recon- 
structed e-machine (solid line). There is generally good agree- 
ment between two, except that the Bragg-like peaks from the 
e-machine are slightly displaced from the experimental spec- 
tra and the maximum from the e-machine overestimates ex- 
periment. The e-machine also has some difficulty reproducing 
the shape of the experimental spectra. The profile 7?.-factor 
between the two spectra is 1Z € m = 29%. 



D. SK229 



FIG. 13: The r = 3 e-machine reconstructed for SK229. 
The large asymptotic state probabilities for the S2 and 5s 
states, as well as the large state-transition probabilities be- 
tween them, show that this is predominately a 2H crystal 
with some faulting. 



Structure Contribution 

2H 82% 

3C 4% 

Deformation fault 9% 

Growth fault 4% 

Other disorder 1% 



Lastly, we examine an as-grown 2H crystal. The 
diffraction spectrum for this crystal is shown in Fig. 1121 
We find the figures-of-merit closest to their theoretical 
values over the interval I e [—0.33, 0.67] with values of 
7 = -0.49 and (3 = 1.00. We find that the smallest-r 
e-machine that gives reasonable agreement between the 
measured and e-machine spectra has a memory length of 
r = 3. The reconstructed e-machine is shown in Fig. 1131 
The large asymptotic state probabilities for the £2 and 
^5 CSs, as well as the large inter-state transition prob- 
abilities between them, indicate this is predominantly a 
2H crystal. More specifically, the probability of seeing 
sequences 1010 and 0101, corresponding to the 2H cycle, 
have a combined total weight of about 82.5%. The re- 
maining probability is distributed among the other four- 
teen length-4 sequences. It is tempting to interpret the 
remaining structure as faults and, indeed, we can. 

Let us treat the transitions s — from 5g and s = 1 
from S\ as though they are missing. These are the 
least probable transitions in the e-machine: P(0,iSe) = 
P(5 6 )P(0|,5 6 ) « 0.004 and P(l,Si) - P(«Si)P(l|,«Si) « 
0.002. Then, in the left half of the e-machine there 
is structure associated with a 2H deformation fault 
[S 7 S 6 S 5 S 3 ] with probability weight 0.040 = P(1011) + 
P(0111) + P(1110) + P(1101). In the right half there like- 
wise is 2H deformation fault structure [<5>o<5i<52<S4] with 
weight 0.049 = ±P(0100)+P(1000)+P(0001)+|P(0010). 
The right portion contains the CSC [52tS4iSi] with weight 
0.036 = iP(0100) + P(1001) + iP(0010), which is asso- 
ciated with growth faults. The CSCs [So] and [Sj], iden- 
tified as 3C structure, have a combined weight of 0.041. 



TABLE II: The Fault Model structural interpretation of the 
reconstructed e-machine of Fig. 1131 This is valid only under 
the assumption of weak faulting. 

Given these observations, the interpretation suggested by 
the e-machine is that SK229 has crystal structures and 
faults in the proportions given in Table ILTl The decom- 
position there is reasonable since there is an underlying 
crystal structure present and the smaller, faulting paths 
are not too large. 

SK229 has not been analyzed quantitatively using the 
FM. By comparing the FWHM of the integer-i to half- 
integer-Z peaks, Sebastian and Krishna^ concluded that 
deformation faulting is the primary vehicle responsible 
for the deviation from crystallinity seen here. We are in 
agreement, except that we also detect small amounts of 
3C crystal structure and some growth faults. 

Figure ITU compares the CFs from experiment and the 
e-machine. The agreement is good, although the recon- 
structed e-machine underestimates somewhat the magni- 
tude of the oscillations in Q s (n). A visual comparison 
of the experimental diffraction spectrum and that gen- 
erated from the e-machine (Fig. I12|) shows that there is 
reasonable agreement. We calculate an 7t!.-factor between 
the two spectra of Tl e M = 29%. 

There are some noticeable differences between the two 
spectra, however. First, we see that Bragg-like peaks in 
the e-machine spectrum are slightly shifted from those in 
the experimental spectrum. Second, the e-machine spec- 
trum overestimates the maximum intensity in each peak. 
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FIG. 14: Comparison of the Q a (n) versus n for experimen- 
tal spectrum SK229 (solid squares) and the e-machine (open 
squares) . 



Third, the peak profiles are qualitatively different. For 
the e-machine, the shoulders of the peaks are broader 
than those of the experiment and the crowns are nar- 
rower. Indeed, the peaks in the experimental spectrum 
appear plateau-like. The sides are very steep and the 
crown is rounded. It is not known if this results from 
instrument resolution or it is an observable effect. 

This example shows that slightly faulted crystals with 
sharp Bragg-like structures can be difficult to analyze. 
Certainly, the basic crystal structure is clear — in this case 
2H, but a very fine experimental /-mesh is needed to map 
each Bragg peak. In contrast, highly disordered spectra 
with significant diffuse scattering are less sensitive to ex- 
perimental details, such as instrument resolution. For 
highly diffuse spectra, where the assumptions underlying 
the FM break down, eMSRboth in practice (extraction of 
correlation information) and in principle (the e-machine 
describes any amount of disorder) reaches its full poten- 
tial. 

From the e-machine we calculate the three length pa- 
rameters to be r = 3 ML, V = 3.5 ML, and A c = 19 ± 2 
ML. The three measures of intrinsic computation are 
found to be = 0.30 bits/ML, C p = 1.8 bits, and 
E = 0.89 bits. 



IV. PHYSICS FROM e-MACHINES 

Now that we have models (e-machines) for the stochas- 
tic processes underlying the observed ML stacking for 
each sample, we can calculate a range of structural, com- 
putational, and physical characteristics that describe the 
stacking patterns and disorder. In Table II I II we list 
the measures of intrinsic computation and characteristic 
lengths calculated for each sample, as well as for three 
crystal structures for comparison. 



A. Characteristic Lengths in Polytypes 

We first note that it was necessary to perform eMSR 
up to r = 3 for each spectrum. This is not surprising 
since the mechanism of deformation faulting is expected 
to be important in ZnS, and the minimum e-machine on 
which deformation faulting structure can be modeled is 
r = 3. This implies, of course, a longer memory length 
that either pure 3C or 2H structure alone requires. The 
generalized periods for the 2H and 3C structures at V = 2 
ML and V = 1 ML, respectively, are also much shorter 
than for the disordered structures which average V = 5.7 
ML. This shows that there is spatial organization over a 
modest range — 6 MLs — for disordered ZnS crystals. This 
pales in comparison to crystalline polytypes with repeat 
distances over 100 MLs, but is still much larger than the 
calculated range of inter-ML interactions of ~ 1 ML. We 
note that many of these long-period crystalline polytypes 
are believed to be associated with giant screw disloca- 
tions that are expressly absent here. For both r and V, 
the disordered structures have values much closer to that 
expected for the 6H structure. In contrast to a perfect 
crystal, the correlation lengths are finite rather than infi- 
nite. Interestingly, the sample that has the most stacking 
disorder (as measured by h^), SK137, also has a compar- 
atively long correlation length. While this was previously 
classified as an as-grown 3C crystal by Sebastian and Kr- 
ishna, 6 we find that it also contains a significant amount 
of stacking sequence associated with the 2H structure. 
Since we cannot assume that any of these structures are 
in equilibrium or the ground state, we can draw no con- 
clusions about the range of inter-ML interactions. 



B. Intrinsic Computation 

Each of the diffraction spectra we analyze also show 
much stacking disorder. Even a spectrum that is quite 
crystalline, like SK134, has a stacking entropy of — 
0.50 bits/ML. Of course, for a crystal the entropy rate is 
zero and for the case of completely disordered stacking 
one would have — 1 bits/ML. If we compare SK134 
and SKI 35, each beginning as a 2H crystal but annealed 
at different temperatures, we see that the latter is slightly 
more disordered, as we expect. The statistical complex- 
ity, a measure of the average history in MLs needed to 
predict the next ML, is also relatively constant at 2.3 bits 
and 2.5 bits, respectively. 

In fact, the measures of intrinsic computation are 
nearly equal, except those for SK229. But SKI 34 and 
SK135 have very different structures. SK134 is largely 2H 
in character while SK135 is largely twinned 3C. Assum- 
ing that they were identical before annealing, this would 
suggest that the disordering process has little effect on 
these measures. This is, of course, tentative, since such 
a conclusion can only be drawn after examining many 
disordered samples. Experimental spectra in the midst 
of the 2H-to-3C transformation would be of significant 
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TABLE III: A comparison of the three characteristic lengths 
and three measures of intrinsic computation that one can 
calculate from knowledge of the e-machine. We calculate 
them for the experimental diffraction spectra, as well as 
for three crystal structures for reference. Recall that A = 
C M - E - rh^ 

interest here. It is possible, though, that SK137 might 
be such a instance. While this is an as-grown twinned 
3C crystal, as noted above this, this crystal also has 
some significant 2H character. Since Sebastian and Kr- 
ishna 6 found that both of these were well described by a 
random distribution of twin faults, they concluded that 
disordered 3C crystals found in the growth furnace re- 
sult from a phase transformation from the 2H structure 
upon cooling the furnace. We find that the two samples 
(SK135 and SK137), while similar, do have qualitative 
differences. We can understand this either as a crystal 
not completely transformed or that the mechanism which 
created SK137 is not simple. We feel that more experi- 
mental data is needed in order to arrive at a more com- 
plete understanding. Since = 0.65 bits/ML for SK137 
and is thus more disordered than either SK134 or SK135, 
the interpretation of this crystal being in the midst of a 
phase transition is a plausible explanation. The most 
striking feature of the measures of intrinsic computation 
is their relative consistency (except for SK229) even while 
the structure of the crystal changes significantly. 



C. Stacking- Fault Energies 



Terms with an odd number of spins do not appear due 
to symmetry considerations. We take Si — ±1 here. 

Engel and NeedsS found that Ji = 0.00187 eV per 
ZnS pair and J 2 = —0.00008 eV per ZnS pair and that 
J3 and K are negligible. Given this let us rewrite Eq. 
in terms of the energy per ZnS pair and take Eq = 0. 
Then the configurational energy is 

E = -Ji(siS i+1 ) - J 2 (s;s i+2 ) , (2) 

where brackets indicate the expectation value over the 
stacking sequence. The expectation values are found di- 
rectly from sequence probabilities, as follows: 

( Sl s l+1 ) = P(11)±P(00) -2P(01) 

(s iS i +2 ) = P(111)±P(101)±P(000)±P(010) 

-2P(110) - 2P(100). (3) 

The configurational energy in terms of meV per Zn-S 
pair is shown in Table HVl for several crystalline structures 
and each of the four disordered samples. The configura- 
tional stacking energies are bounded above and below by 
the 2H and 3C structures with relative configurational 
energies of 1.95 meV/ZnS-pair and -1.79 meV/ZnS-pair, 
respectively. For SKI 34, the annealing process has in- 
troduced faults and has lowered the average stacking en- 
ergy from the original 2H structure to 1.13 meV/ZnS- 
pair, while increasing the stacking entropy. If we assume 
that SK137 is a partially transformed 2H-to-3C crystal 
(though mostly 3C), then we see that the crystal expe- 
riences further disordering and the stacking energy falls 
to -0.57 meV/ZnS-pair. SK135 shows the most advanced 
transformation with the 2H structure almost completely 
eliminated and stacking energy not too far from the ideal 
minimum at -1.02 meV/ZnS-pair. The stacking entropy 
begins to fall, however, as the transformation nears a dis- 
ordered 3C crystal. Being only a slightly disordered 2H 
crystal, SK229 shows the highest stacking energy of 1.56 
meV/ZnS-pair. As we might expect from the relative 
magnitudes of J\ and J 2 , the contribution from the J\ 
term completely dominates the energy. 



One physical quantity amenable to calculation from 
the e-machine is the difference in configurational energies 
of the particular polytypes. Numerical calculations find 
that the configurational energy depends only the nearest 
and the next-nearest neighbors in the stacking arrange- 
ment. Engel and Needs^ have done a first-principles 
pseudopotential calculation of the total energy of five ZnS 
polytypes, from which they determined the strength of 
the interactions up to the third nearest layer. The most 
general expression possible for inter-ML interactions up 
the third nearest neighbors is given byii 

E = Eq - Ji 2J SjSj+l - J 2 SiSi+2 

i i 

-J3 / ] SjSj+3 - K)^SiSi + iSi +2 Si + 3 . (1) 

i i 



D. Hexagonality 

The degree of birefringence of ZnS crystals is known to 
depend only a single structural parameter — the hexago- 
nality a/t^ 6 - This parameter is defined as that fraction 
of MLs which are hexagonally related to their neighbors. 
That is, a/i is defined as the frequency of occurrence of 
sequences ABA and BAB and their cyclic permutations. 
In terms of the Hagg notation, these are simply P(01) 
and P(10), respectively. Since P(01) = P(10),i we have 

a h = 2P(01) . (4) 

Sequence probabilities are directly calculable from the e- 
machine, so that the hexagonality of a disordered crystal 
can be easily found. 
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System 




{SiS i+2 ) 


E ah History 


2H 


-1.00 


1.00 


1.95 1.00 PC 


3C 


1.00 


1.00 


-1.79 0.00 PC 


6H 


0.33 


-0.33 


-0.65 0.33 PC 


SK134 


-0.58 


0.63 


1.13 0.80 D 2H, 300 C for lh 


SK135 


0.56 


0.45 


-1.02 0.24 D 2H, 500 C for lh 


SK137 


0.32 


0.45 


-0.57 0.34 AG 3C 


SK229 


-0.80 


0.86 


1.56 0.90 AG 2H 



TABLE IV: Relative configurational energies E of experimen- 
tal polytypes and several pure-crystalline polytypes. The last 
column gives the history of each sample, where PC stands 
for perfect crystal, AG as-grown, and D disordered. We use 
the energy coupling constants, Ji and J2, calculated by En- 
gels and Needs along with the reconstructed e-machine for the 
disordered processes to find the configurational energy of the 
disordered structures via Eqs. and Q. 

In Table IIVI we show the hexagonality calculated for 
all of the spectra, as well as for several crystal structures 
for comparison. 



V. DISCUSSION 

We have successfully applied computational mechan- 
ics to the discovery and description of stacking order in 
single crystals of polytypic ZnS. In doing so, we recon- 
structed from experimental diffraction spectra the mini- 
mal, optimal, and unique description of the stacking pro- 
cess as embodied in the e-machine. In contrast to pre- 
vious analyses^ we used all of the information in the 
diffraction spectra, both in the Bragg peaks and in the 
diffuse scattering between them. We imposed no restric- 
tions on the kind of structures to be found, save that they 
be representable by e-machines. Further, the computa- 
tional mechanics approach was not limited to the case of 
weak faulting, but can be used to treat even highly disor- 
dered samples. Additionally, the e-machine can naturally 
accommodate more than one parent crystalline structure 
as seen in SK134. 

For two of the spectra, a sensible decomposition^ of 
the e-machine into crystal and faulting structure was pos- 
sible, allowing a direct comparison between the compu- 
tational mechanics approach and the FM. The e-machine 
detected structures not previously found by the FM. For 
example, in SK134 we found that not only was struc- 
ture associated with deformation faulting important, but 
there was also structure related to growth faults and 
layer-displacement faulting. We even found nascent se- 
quences leading to the 3C structure. For the other two 
cases, while no FM-like decomposition of the e-machine 
was proposed, we still found significant structure as em- 
bodied in the e-machine. From the e-machine, physical 
insight into the structure of the stacking was possible. 
For example, in the r = 3 reconstructed e-machine for 
SK137 we could eliminate 6H structure based on the ab- 



sence of a transition between CSs. We also found that 
2H structure was present. Even when no sensible de- 
composition into a simple pure-crystal and weak-faulting 
structure is possible, the e-machine still directly provides 
sequence frequencies, which can be used to build physical 
insight into the stacking structure. 

From a knowledge of the e-machine, it is possible to cal- 
culate a number of physical characteristics. In Table HVl 
we tabulated the stacking entropy per ML for each spec- 
trum. Given the coupling parameters between MLs we 
calculated the average stacking-fault energy for the sam- 
ples, as shown in Table IIVI We were also able to find 
the hexagonality for the disordered crystals. Knowing 
the e-machine allowed us to find various characteristic 
lengths associated with each disordered crystal, such as 
the memory length and the generalized period. We be- 
lieve that additional physics can be calculated from the 
e-machine. 

We also calculated measures of intrinsic computation 
from the e-machine. We found that the minimum mem- 
ory length for all samples was r = 3, which is in excess of 
the calculated inter-ML interaction range of ~ 1 ML for 
ZnS. We further found that the statistical complexity, a 
measure of the amount of information in the e-machine, 
was also much larger than that of either the pure 2H or 
3C structures. Also, the range over which structures are 
found in the disordered samples is about 6 ML. 

Characterizing solid-state transitions in polytypic ma- 
terials is of considerable interest. Let us review what an 
e-machine does and does not imply. Most simply put, 
the e-machine is the answer to the question, What is 
the minimal, optimal, and unique description of the one- 
dimensional stacking structure of the sample? Any phys- 
ical parameters that depend on this description are in 
principle calculable from e-machine. The e-machine does 
not answer the question, How did the crystal come to 
be stacked in this way? To determine this, one must 
augment structural knowledge (as embodied in the e- 
machine) with additional information or assumptions. 
Such information can come in the form of a time series 
of structures obtained either from a series of numerical 
simulations or experiments or, in the theoretical domain, 
perhaps from assumptions about weak faulting. 

Since we have discovered structures in polytypic ZnS 
that were undetected before, we feel that the mecha- 
nism of faulting — previously attributed to deformation 
faulting — deserves re-examination. We have provided a 
firm theoretical foundation for the discovery and descrip- 
tion of disordered stacking sequences in polytypes and be- 
lieve that additional experimental studies are warranted. 
Additionally, computer simulations of solid-state trans- 
formations in polytypes with proposed faulting mecha- 
nisms, accompanied by the concomitant reconstruction 
of the e-machine directly from the sequence of MLs from 
the simulation, should provide a powerful method to un- 
derstand the gross features of the transformations. 

Before such studies can be definitive, a quantitative 
understanding of the effects of experimental error on the 
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reconstructed e-machine is needed. With the exception 
of our introduction of the figures-of-merit, 7 and /3, we 
have not addressed this important issue here. We note 
that the original experimental data was not reported with 
error bars and this means that comparison with the de- 
sired e-machine error analysis would have not been pos- 
sible. Additionally, the necessity to digitize the data un- 
doubtedly introduced errors. It is therefore difficult to 
assess the amount of error in each spectrum. Our intu- 
ition tells us that error in the diffraction spectrum will 
likely lead to suppression of the more delicate structures 
on an e-machine. Therefore, one should expect to find 
less structure and more randomness. 

We mention that the application of computational me- 
chanics to the description of one-dimensional sequences 
is the most general approach possible to this problem. 
Thus, its application here to polytypism represents the 
end point in the evolution of models to describe the dis- 
ordered sequences seen these substances. Any alternate 
description can be expressed as an equivalent e-machine 
and none can be more general, since in the language of 
statistics an e-machine is the minimal sufficient statistic 
for the underlying process. It may be possible to find 
specialized algorithms that are more sensitive or more 
efficient in determining an e-machine than the one we in- 
troduced in e-machine spectral reconstruction. However, 
the answer, in its most general form, will be expressible 
as an e-machine. 



We also note that our work here represents a solution 
to a significant theoretical problem — How does one ex- 
tract structural information from a power spectrum? Our 
application has been to polytypism, but the principles 
underlying our solution may be applied quite generally 
to domains in which spectral information is available. 

Future directions for this work include an application 
of eMSR to other polytypes, as well as to substantially 
more complex materials. The extension of these ideas to 
the more common cases of disorder in two and three di- 
mensions is also desirable. The development of computa- 
tional mechanics in higher dimensions would significantly 
aid in the classification and understanding of disorder in 
many physical systems. 
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